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A HEAVY ION BEAM PROBE SYSTEM FOR INVESTIGATION OF A 
MODIFIED PENNING DISCHARGE 
by George X. Kambic * and Walter M. Krawczonek 
Lewis Research Center 

SUMMARY 

An ion beam probe diagnostic system can measure time- and space-resolved pro- 
files of plasma space potential and electron density. In combination with a computer 
iterative technique, the ion beam probe can determine both the space potential profile 
in plasmas containing strong electric fields and potentials comparable in magnitude to 
the energy of the probing ion beam. During ion beam probing of a modified Penning dis- 
charge, several groups of secondary ions have been observed coming from the plasma 
with a fixed primary-beam energy and momentum. The energies of these ions were 
within 10 percent of the values predicted by a computer-generated model of the potential 
profile in the plasma. The mechanical and electronic components of the system are 
described, with particular emphasis on those features required to probe plasma poten- 
tials comparable in magnitude to the ion beam energy. 

INTRODUCTION 

The heavy ion beam probe (refs. 1 and 2) is a nonperturbing plasma diagnostic. It 
has the unique ability to measure plasma space potential directly. For this reason, it 
is important to develop this diagnostic for the investigation of laboratory- scale and 
fusion- related plasmas containing strong electric fields (ref. 3). It can make time- and 
space -resolved measurements of electron density, electron temperature, plasma po- 
tential, and axial current density (refs. 4 to 6). Interest in this diagnostic method at 
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Lewis stems from the existence of plasmas containing strong electric fields in their 
interiors. These devices include the modified Penning discharge (refs. 7 and 8) and the 
bumpy torus (refs. 9 and 10). This report describes the mechanical and electronic sys- 
tems required to implement ion beam probing on the modified Penning discharge. 

The probe built at Lewis is based on the ion beam probing technique developed by 
Hickok and Jobes (refs. 4 and 5). Figure 1 illustrates the principle of the ion beam 
probe. A monoenergetic ion beam (f 1 ") is injected into a plasma confined by an axial 
magnetic field. A fraction of this primary beam, called the secondary beam, is doubly 



Figure 1. - Principle of ion beam probing spatial resolution and space potential measurement. 


ionized and magnetically separated from the primary beam and is detected outside the 
plasma volume. Secondary ions are generated along the length of the primary-beam 
path through the plasma, but a given detector sees only those generated in a small vol- 
ume of the plasma. By measuring the properties of the secondary beam, including en- 
ergy, current, and momentum, the previously listed plasma characteristics can be deter- 
mined. In this report, we concentrate on measuring the plasma potential profile. The 
potential is determined by calculating the energy of an ion as it passes through the 
plasma. The ion has an initial energy Eq at s = 0, where s is the length along the 
beam path: 


e(s )| b=0 = e(°) = Eo 


(i) 


2 



As the ion penetrates the plasma, it changes total energy as determined by the space po- 
tential cp : 

b 


ECSj) = E q - e<p s (2) 

If the ion is then ionized to a +2 state, the final secondary- ion energy is equal to the 
initial primary- ion energy plus the plasma space potential, as given in equation (3). 


E(s 3 ) - E final 

(3a) 

E(s 3 ) = E q - ecp s + 2ecp s 

(3b) 

E(s 3 ) = Eq + e<p g 

(3c) 


-2 -3 

In previous ion beam probing, plasma potentials have been about 10 to 10 times the 
primary-ion-beam energy, ecp « E^ In the modified Penning discharge plasma, the 
potential can be many kilovolts (ref. 8). With the Lewis probe, we measure potentials 
that are large fractions, 0. 25 to 0. 50, of the primary-ion-beam energy. 

The Lewis ion beam probe uses a monoenergetic beam of thallium ions (Tl + ). The 
beam is injected into the midplane of the modified Penning discharge, which is sustained 
by a pair of Penning anode rings. The maximum magnetic field used was 0. 46 tesla. 

For a thallium ion (mass, 204. 4 atomic mass units) of 20-keV energy, the ion trajec- 
tories through the magnetic and electric fields of the plasma cannot be approximated by 
a simple trajectory appropriate to motion in the magnetic field alone. If the maximum 
potential in the plasma is of the order of the anode voltage, the secondary-ion energy 
can change by 50 percent. Since the ion trajectories are no longer uniquely known, 
spatial resolution is lost. 

This problem was dealt with by using a computer program to calculate a self- 
consistent model of a plasma potential profile. In its simplest form, it solves the equa- 
tions of motion of a particle of mass m and charge q in electric and magnetic fields: 

mx = q(E + x x B) (4) 

The magnetic field is known everywhere and the electric field is obtained by differentia- 
tion of the potential profile. This calculation is described in more detail elsewhere 
(ref. 11). The shape of the potential is varied in an iterative manner until measured and 
calculated trajectories are in agreement. 
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ION BEAM PROBE APPARATUS 


The Lewis ion beam probe uses a thallium ion beam with energies of 15 to 40 keV. 
The modified Penning discharge facility is shown in figure 2. The beam is injected into 



Figure 2. - NASA Lewis modified Penning discharge and ion beam probe facility. 


the midplane of the mirror magnetic field of the modified Penning discharge. The large 
vacuum tank contains the Penning anode rings and the magnetic field coils. A cutaway 
drawing of the plasma region is shown in figure 3. The Penning anode rings are shown, 
along with the grounded cylindrical wire mesh that surrounds the plasma and symmet- 
rizes the radial potential in the midplane of the system. Figure 4 shows the screen, the 
anode support shaft, and one of the two superconducting coils. The plasma is confined 
to the region inside the anode rings. The ion beam passes between the double anode 
rings perpendicular to the magnetic field. The primary beam is detected with either the 
electrostatic analyzer or the primary detector, and the secondary beam is detected with 
the electrostatic analyzer. 

Figure 5 is a schematic of the various components of the plasma and probe system; 
it shows the ion gun, the anode ring, the grounded screen, the electrostatic analyzer, 
and the primary detector. The primary detector is a set of 33 flat Langmuir probes 
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Figure 3. - Cutaway drawing of apparatus. 


Figure 4. - Magnet Dewar, anode ring, and grounded screen. 
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placed on the tank wall in the path of the primary beam. The overall deflection angle of 

the beam going into the electrostatic analyzer is 40°. The maximum midplane magnetic 

field of which the coil system is capable is 0. 8 tesla, but studies were done at a con- 

2 2 

stant 0. 46 tesla. Since the scaling of the ion beam can be expressed as E^rm °c B r , 
the beam probe characteristics of other systems can be scaled from the characteristics 
of this system. 

The ion gun used is shown in figure 6(a). It is basically the same gun used by 
Hickok and Jobes (ref. 5). It uses Pierce extraction electrodes near the ion source 
(ref. 12). It is operated by means of a resistive divider string and keeps the beam full 
width at half maximum to within 0. 5 centimeter over the beam energy range used. The 
optimum voltage ratios for the electrodes of the gun were found by trial and error. The 
gun is dependable and convenient in operation. 

The thallium ion source is a thallium zeolite packed into the top of a 0. 64- 
centimeter-diameter ceramic tube, shown in figure 6(b). The tube has an 0. 5- 
centimeter-diameter insert that is inset 0. 64 centimeter. A 0. 025-centimeter-diameter 
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(b) Details of thallium ion source. 
Figure 6 . - Electrostatic ion gun. 
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tungsten filament is placed in the well at the top of the ceramic tube for thermionic 
heating of the zeolite. Tantalum leads 0. 05 centimeter in diameter are welded to the 
tungsten filament. The sources were run for 40 to 50 hours at constant beam current. 
The zeolite material can be doped with other materials such as sodium or lithium. 

The electrostatic energy analyzer is a parallel-plate, energy selection system that 
operates around a mean entrance angle of 45° and a geometric gain of 1 (refs. 13 to 15). 
At this angle, the voltage across the analyzer in volts is numerically equal to the energy 
of a doubly charged secondary ion in electron volts when the height of the analyzer is 
equal to the separation of the ion entrance and exit slits. A diagram of the analyzer is 
shown in figure 7. This analyzer is sensitive to small changes in the beam energy. It 
was useful in detecting an 8 -volt oscillation in the ion gun power supply when the 
primary-ion energy was 6 keV. At this time the analyzer was detecting primary 
ions so the voltage across the analyzer V ga was 12 keV. Although the system can 
measure small changes in beam energy, we used it to measure plasma potentials that 
were many kilovolts. 

The location of the electrostatic analyzer shown in figure 5 was chosen on the basis 
of calculated ion trajectories. These trajectories took into account a potential profile 
based on a model for the modified Penning discharge developed by Roth (ref. 8). The 
split-plate detector of the +2 ions can operate as an infinitely thin line detector when 
used in conjunction with a feedback amplifier. The feedback circuit is part of the elec- 
tronic system for the probe. The associated electronic system is shown in a block dia- 


Analyzer 

structure 



Figure 7. - Electrostatic energy analyzer and feedback loop. 
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gram in figure 8(a) and in schematic form in figure 8(b). One innovation in this system 
is the use of separate high-voltage power supplies to operate the ion gun and the analy- 
zer, a result of the high space potentials existing in the plasma. 

The electronic system comprises (1) the scan control and variable-time-base gener- 
ator, (2) the initial value sets (for the energies and sweep angles) and scaling (for the 
ion mass and magnetic field), (3) the feedback circuitry and beam detection, and (4) data 
acquisition. 

One function of the electronic system is to control the primary-ion-beam energy and 
injection angle into the plasma. The initial ion gun energy can be set, and the magnitude 
of the magnetic field can be scaled. The secondary beam is detected in the split-plate 
detector and amplified in a fast amplifier. The feedback amplifier for the detection sys- 
tem serves to focus the secondary ions onto the detector by varying the voltage across 
the analyzer. In a plasma where the space potential is small compared with the ion en- 
ergy, this double power supply arrangement is unnecessary and the plasma potential can 
be measured in the manner described in references 1 and 2. 

The scan control and variable-time-base generator starts and stops beam sweeps 
through the plasma, sets the number of ionization points that are to be measured, and 
controls the length of time the primary beam stays at each ionization point. The timing 



(a) Block diagram. 

Figure 8. - Electronic system associated with feedback circuit 
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Figure & - Continui 
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(b) Continued. 
Figure a - Continued. 
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circuitry provides signals to a deflection step control that causes the energy across the 
ion gun resistive divider to change, thus changing the energy of the beam injected into 
the plasma. It also changes the voltage across the deflection plates of the ion gun, thus 
varying the injection angle of the primary beam into the plasma. Since varying the 
primary-beam energy and injection angle moves the ionization point around in the plasma, 
an array of points inside the plasma volume is sampled, each of which is an ionization 
point in the ideal case. The number of ionization points examined can be controlled by 
appropriate scaling of the scan control and by adjusting the increments in both energy 
and angle. The scan control also provides a signal to the power supply that drives the 
electrostatic analyzer. 

The secondary ion beam results from the interaction of the plasma with the primary 
beam. The secondary beam enters the analyzer and is deflected into the split-plate de- 
tector, where the total current is measured. A detector amplifier measures the current 
and provides a signal to a feedback circuit that drives the analyzer power supply in order 
to keep the secondary beam focused on the split-plate detector. The voltage across the 
analyzer is varied until the current to each plate of the detector is equalized. Since the 
feedback signal to the analyzer is proportional to the difference in the energy of the pri- 
mary ion beam as compared with the energy of the secondary beam, this feedback signal 
is a measure of the plasma potential. The total secondary-beam current, which is pro- 
portional to the electron density, can also be recorded. The feedback signal or the total 
secondary current signal is sent through a voltage-to-frequency converter and is then 
stored in an appropriate channel in a multichannel analyzer. The magnitude of each sig- 
nal can be read out on the z-axis of an oscilloscope (intensity), while the x- and y-axes 
can be related to x- and y-positions inside the plasma. The data can also be printed or 
punched out for further computer analysis. 

This control system is basic to any ion beam probing system that is used for time- 
or space-resolved plasma measurements. This particular system has not been used as 
a complete driving and feedback control system because of time limitation and an inade- 
quate secondary-beam signal. The secondary signals are always less than 1 nanoampere 
and have usually been less than 100 picoamperes. These signals can be detected but are 
too small to be used in the present secondary-beam detector amplifier circuit. The 
electron density implied by the magnitude of these signals is approximately 10*° to 10 11 
per cubic centimeter. There is usually a unique secondary signal associated with each 
set of primary-beam initial conditions. However, for this plasma, we have found that 
more than one secondary beam can be observed in the electrostatic analyzer for one set 
of primary-beam conditions. 

The expected current can be calculated from the relation between the primary and 
secondary ion beam currents: 
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(5) 


— = 2N i2TLl 

hi % 

where 1^ and are the primary- and secondary-beam currents, respectively, 

is the velocity of the beam ions, N g is the electron density, and (ov)/^ is the effective 

cross section (ref. 5). The term ( ov) is obtained by averaging the cross section for 

ionization over an assumed Maxwellian electron distribution, and l is a resolution 

length of the electrostatic analyzer along the beam. For this example, a good estimate 

for 2 is 1 centimeter, which is obtained by varying the ionization point of the secondary 

beam and observing the effect at the electrostatic analyzer. The secondary currents 

-5-3 

would then be estimated as 10 to 10 microampere (using 1^ = 0. 5 ji A, Eq = 20 keV, 
(crv}/v^ ~ 10” 14 cm^ (ref. 5), at an estimated electron temperature of 100 eV). 


COMPUTER CALCULATIONS 

The modified Penning discharge has strong electric fields in the interior of the 
plasma. The computation must correctly locate the ionization point at which the ob- 
served secondary ions are being emitted. The location of the ionization point must be 
compatible with the ion beam initial conditions, the observed secondary ion energies, 
and the plasma space potential profile. 

The calculation of a particle trajectory is a straightforward problem, which is found 
by solving equation (4). This is the equation of motion of a particle of mass m and 
charge q in electric and magnetic fields, E and B, respectively. The magnetic field 
is specified everywhere along the orbit. This is true for low-/3 plasmas, such as the 
modified Penning discharge where p is less than 10~ 4 (p is the ratio of plasma pres- 
sure to magnetic field pressure (ref. 16)). Also, the current density in the plasma must 
be low. In the modified Penning discharge, the total current is less than 1 ampere. A 
current along the magnetic field axis would create a poloidal field that would bend the 
beam out of the two-dimensional plane used to model this plasma. It is possible to 
measure the current density by using this deflection (refs. 17 and 18). 

The computer program input consists of the primary-beam initial conditions and the 
model potential profile parameters. For the primary-beam trajectory calculation alone, 
this is enough information. The calculation of the best-fitting parameters can then pro- 
ceed as described in reference 11. Calculating the secondary orbits is no more complex 
than calculating the primary orbits. The program picks points along the primary path 
where the charge of the ion is set equal to 2e. The actual ionization point in the plasma 
is unknown, but there is a large change in the secondary ion energy after it leaves the 
plasma. The program uses this to locate ionization points. The program searches for 
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the largest change in the secondary ion energy. Using the model potential, it attempts 
to match this energy with a particular measured secondary ion energy. It does this by 
continuing to vary the model parameters in an iterative loop. 

The anode rings and grounded screen in the modified Penning discharge are concen- 
tric. It is convenient to assume circular symmetry for the model and to ignore effects 
due to the anode support shaft. With the symmetry assumption we can numerically solve 
for the vacuum field by using a program that solves Laplace's equation (ref. 19). This 
vacuum field is shown in figure 9 along with a plot of the assumed potential (ref. 8). It 
provides a first step to the model profile. When the plasma is present, the anode is 



Figure 9. - Initial potential profile and vacuum field. 


shielded within a Debye length, and the potential is expected to be less than the vacuum 

field inside the anode ring. Previous research on the discharge indicates a model for 

the shape of the potential inside the anode (ref. 8). The particular potential obtained is 
o 

a function of r and, if the potential is assumed never to be less than zero, it can be 
written in the form 

<p = CjT 2 + C 2 (6) 

o 

where (p = 0 when r = -Cq/C^. A potential of this type is assumed for the region in- 
side the anode ring. We model the potential outside the anode by the vacuum field and 
inside by the function in equation (6). Dow (ref. 20) has suggested that there may be 
potentials greater than zero on the axis of a Penning discharge. The model can allow 
for this by scaling the potential in equation (6). Finally, a flat region is inserted near 
the peak potential. The resulting model is shown in figure 10. The five parameters are 
^max’ ^in’ R i’ R o’ and R r' In the re S ion R o- — r — the potential is assumed to 

o 
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Figure 10. - Assumed model of potential profile. 

be a linear function of r. As can be seen from” the vacuum field solution in figure 9, 
this is a good approximation. The potential assumes the value <P max in the region 
R q > r > IT. The potential is described by equation (6), and this equation is scaled 
depending on the value of <p- n in the region > r > R r> The potential would be di- 
vided into as many regions as needed to improve the fit of the final potential. The com- 
putations are one dimensional. 


EXPERIMENTAL PROCEDURES 

Experimental results obtained in this device have been reported elsewhere (ref. 11). 
A summary of the experimental methods and procedures is given here. 

The primary beam was first examined with no plasma present. The results are 
given in table I. One important point of this table is that, for the single magnetic field 
strength used, a beam energy of about 6 keV was sufficient to deflect the primary beam 
into the electrostatic analyzer. It was found that a systematic error of about 1.50 de- 
grees exists in the computation of the orbit. Estimates of the error in the location of 
the beam emission point caused the overall estimate of angular error to be about ±2 de- 
grees. For these initial measurements, this was considered acceptable. Table EE gives 
data on the primary-beam trajectories for various discharge conditions. Listed are 
the primary-beam energy, injection angle, and trajectory end point and the plasma dis- 
charge conditions. Beam angles varied from 43. 3° to 35. 3°. The plasma is character- 
ized by the Penning anode voltage and the background deuterium pressure, which is 
measured with an ionization gage. The voltage ranges from 4. 5 to 17 kilovolts, and the 
background neutral density from 2. 8x10^ to 1. 6x10^ neutral particles per cubic centi- 
meter. The best parameters, in terms of obtaining usable experimental data, were an 
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TABLE I. - TYPICAL MEASURED AND CALCULATED 


PRIMARY-ION-BEAM TRAJECTORIES WITH 


DISCHARGE TURNED OFF 


Primary beam 

End point 

Type of result, 
measured (M) 
or calculated (C) 

Energy, 

keV 

Injection 

angle, 

deg 

Radius, 

r, 

cm 

Angle, 

9 + 0. 5°, 
deg 

6. 15 

39. 8 

92. 5 

1.6 

M 

6. 1 

39. 5 



. 8 

C 

6. 1 

40 



1.0 



6.2 

39. 5 



1. 1 



6.2 

40 


r 

1.3 



20 

41 

42 

7 

25.4 

M _ 

20 

41 



26.6 

C 


22 

42.2 



25.4 

M 

22 

42 



24. 5 

C 


22 

42. 5 



24.6 

c 


17. 8 

42. 5 



25.4 

M 

17. 8 

42.5 



22.2 

C 


31.8 

39. 5 



22.4 

M 

31.8 

39. 5 



20.3 

C 


18 

43. 1 



25.4 

M 

18 

43. 1 



22.4 

C 



anode voltage V 0 of approximately 10 kilovolts and a background neutral density of 
11 a 11 

4. 5x10 to 5. 5x10 neutral particles per cubic centimeter. Problems encountered 

under other conditions are discussed in the next section. 

Table II points up some unique problems of probing a high-voltage plasma discharge 

with an ion beam. For low discharge voltages (~5 kV) the trajectories are governed 

primarily by the magnetic field. Examples of this are beam energies E. of 19. 8 and 

10 8 1 

30. 6 keV in table II. For plasma operating conditions of 1. 7x10 /cm < 

N q < 5. SxlO^/cm 3 (0. 5x10"® torr < P < 16x10"® torr) and anode voltages above 5 kilo- 
volts, the effects of the electric field became apparent. For energies between 11.5 and 
22. 0 keV, the beam was deflected through an angle of 39. 3°, in order to be observed in 
the electrostatic analyzer. A beam deflection much less than 39. 3° could be seen on 
the primary detector. A plot of a trajectory of these beams shows that the electrostatic 
field was changing the curvature of the orbit due to the magnetic field, as illustrated in 
figure 11. 

The experimental results given in table II were examined by using an orbit compu- 
tation program with an assumed potential profile. An initial estimate of the radial 
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Figure 11. - Calculated trajectories for model potentials. 






TABLE II. - EXPERIMENTAL PRIMARY-ION-BEAM INITIAL 


CONDITIONS AND END POINTS FOR VARIOUS 


PLASMA OPERATING CONDITIONS 


Primary beam 

End point 

Plasma 

Energy, 

Injection 

Radius, 

Angle, 

Anode 


Pressure, 

P 

keV 

angle, 



9 + 0. 5°, 

voltage, 

Neutral 

torr 


deg 

cm 

deg 

v a> 

particles 








kV 

per cm^ 



17. 8 

36.9 

42 

7 

22.4 

10.4 

5. 5x10 11 

15.9X10' 6 

17. 8 

43. 1 

92. 5 

1.6 

10.4 

5. 5 

15.9 

22.0 

38.5 

42.7 

25. 4 ' 

13.2 

2 

8xl0 10 


8 

22.0 

42.9 

42.5 

1.6 

13.2 

2 

8 


8 

17.6 

36.9 

42.7 

22.4 

10 

5 

5xl0 n 

15.9 

17.6 

43 

92. 5 

1.6 

10 

5.5 

15.9 

17.3 

38.3 

42 

7 

22.4 

10. 1 

1 

7x10™ 


5 

17.3 

43. 1 

92.5 

1. 6 

10. 1 

1 

7 


5 

11. 5 

36 

42 

7 

22.4 

6 

5. 5X10 11 

15.9 

11.5 

42.6 

92 

5 

1.6 

6 

5 

5 

15.9 

31.7 

39.3 

42.7 

22.4 

10 

5. 5 

15.9 

18 

37. 7 

42 

7 

25.4 

9. 4 

4.6 

13.25 

24. 5 

36.9 

42.7 

25.4 

13.2 

2. 8*10^ 



23.6 

42.9 

92. 5 

1.6 

17 

4. 6xl0 ]1 

13.25 

23. 1 

37.6 

42.7 

20.9 ‘ 

7. 5 





12. 5 

34.7 

42.7 

25.4 

5. 8 





16 

47.86 

92. 5 

1.6 

12 





18 

35.4 

42 

7 

20.9 

10 





18 

35.2 

42. 7 

22.4 

10 





19.4 

41.7 

92. 5 

1.6 

9.8 





19.8 

44.3 

42.7 

22.4 

4.5 





19. 1 

36.8 

42 

1 

22.4 

10 

4 

7 

13 

5 

20.3 

41 

92.5 

1.6 

11.2 

4.6 

13.25 

20.3 

36.7 

42 

7 

25.4 

11.2 





30 

38.6 



20.9 

11.2 





30.6 

42.7 



25.4 

5.3 





37.2 

38.9 



22.4 

5.4 

9. 1x10™ 

2 

.65 

31.6 

41. 1 



26.8 

5 

1.8 

53 



electrostatic potential is used in the program. The parameters shown in figure 10 
divide the radial potential profile into four regions. The vacuum field is a computed 
solution of Laplace's equation in cylindrical symmetry, which is assumed for the 

modified Penning discharge. Inside the radius R., the model potential is proportional 

2 1 
to r plus a constant, which gives an electric field proportional to r in this region 

(ref. 8). 

To optimize the model parameters, a minimization calculation is carried out, with 
the quantity M defined as 


M = £ «> m l - 9 </ n=l,2,... (7) 

i= 1 

where the angles in this expression are the coordinates of the trajectory at the detector 
and © mi and are from the measured and calculated beams, respectively. The 
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calculated angles are a function of the five parameters (R., R q , R r , and <P in )- 

The minimization involved a choice of an initial estimate and the iterative variation of 
the five parameters. The order of variation was <P max > R r , <P in , R p and R 0 - The 
initial estimate is shown in figure 11 and is consistent both with the vacuum field solu- 
tion for r > R = R„ and with the discussion of reference 8 for R„ ^ r < R. = R„. 

The primary-beam initial conditions and trajectory are taken from the first two columns 
in table III. The initial estimate for the two cases examined in detail is specified by 
R. = R q = R a , R r / R a = 0- 87, <p in = 0, and <P max — 0- 88 V a - The iteration processes 
used are described in reference 11. The results of the calculation are given in figure 12 



Figure 12. - Potential profiles for cases 1 and 2 (table III), obtained 
by minimization of the parameter M (eq. (7)). 


as radial potential profiles. The characteristics of the potentials are (1) V max is close 
to V a and (2) there are nonzero potentials on the discharge axis. The result that 
R. > R is corroborated by two pieces of evidence: (1) the location of R./R > 1 and 
V R a > R j/ R a minimizes M, and (2) this potential predicts multiple secondary en- 
ergies that are observed in the experiment. Within experimental error, the potentials 
in the modified Penning discharge are a maximum close to the anode voltage and non- 
zero on the axis. 

Unless the ion orbits cross the discharge axis, we cannot have direct knowledge of 
the region inside the distance of closest approach to the plasma center. The profiles 
in figure 12 were used to predict secondary-ion trajectories from points along the 
primary-ion beam. The result was that under some plasma conditions (e.g. , those in 
case 1 in table HI) multiple secondary-ion beams are created and can be detected in 
the electrostatic analyzer. The potential from case 1 of figure 12 was used in this cal- 
culation. The energy of the primary ions is specified. If the potential is that given in 
figure 12, secondaries emitted from the primary beam at three points in the plasma 
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TABLE IH. - COMPARISON OF CALCULATED AND MEASURED 


END POINTS FOR TWO DISCHARGE CONDITIONS 


Case 

Primary beam 

End point 

Type of result, 
measured (M) 
or calculated (C) 

Plasma 

Energy, 

keV 

Injection 

angle, 

deg 

Radius, 

r, 

cm 

Angle, 

s, 

deg 

Anode 

voltage, 

v a> 

kV 

Pressure, P 

N/m 2 

torr 

1 

17. 8 

36.9 

42.7 

22. 4 

M 

10.4 

5 5X10 10 

15. 9X10' 6 




42.7 

22.2 

C 






43. 1 

92. 5 

1.6 

M 

10.4 

55 

15.9 




92. 5 

8.9 

C 




2 

22.0 

38. 5 

42. 7 

25.4 

M 

13.2 

2. 8X10 10 

0. 8xl0’ 6 




42. 7 

25. 1 

C 






42.9 

92. 5 

1.6 

M 

13.2 

2.8 

o;.8 




92. 5 

3.9 

C 





would enter the electrostatic analyzer. The energies and the emission points are speci- 
fied in table IV and are also shown in figure 13. Such peaks were observed experimen- 
tally for the primary beam and plasma conditions listed in figure 13. Table IV and fig- 
ure 13 show that we have established a possible potential for the plasma. The primary 
beam is observed in the electrostatic analyzer as predicted by the computer. Addi- 
tional adjustment of the model potential profile is required to minimize the difference 
10 percent) between calculated and measured orbits for the secondaries. 


TABLE IV. - CALCULATED SECONDARY- 
ION-BEAM TRAJECTORIES AND 
ENERGIES FROM CASE 1 OF 
TABLE m AND COMPARISON 


WITH EXPERIMENT 


Calculated 

Calculated 

Measured 

energies, 

ionization point, 

energies, 

keV 

(x, y), 

keV 

■ ■ - ■ ■ 

cm 



25 

-9. 2, -6.0 

23 

27 

-8.2, -5.0 

26. 5 

22 

5.2, 3.0 

19.8 
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y, cm 


Figure 13. - Predicted secondary-beam paths for case 1 (ta. 



DEVELOPMENTAL PROBLEMS 

High-voltage breakdown is a problem when the plasma being probed emits energetic 
particles and radiation. The modified Penning discharge produces hot ions, ultraviolet 
rays, and soft X-rays. The hot ions produce noise in the electrostatic analyzer. De- 
pending on the discharge background pressure and the anode voltage, a signal may be 
received in the analyzer with the ion gun off. These ions originate in the plasma. Since 
they are most probably D + or Dg + , they are observed in the analyzer at voltages nu- 
merically equal in volts to twice the ion energy in electron volts. Hence, a 10-keV ion 
would be observed at 20 kilovolts. These ion energies are equal to or slightly less than 
the anode voltage. The intercepted current in the analyzer from these plasma ions is 
large enough to mask the secondary current from the ion probing beam. A simple sec- 
tor magnet or a second set of electrostatic deflector plates was used to filter out these 
plasma ions. 

The flow of background plasma into the electrostatic analyzer and the ion gun must 
be kept to a minimum, for example, by differentially pumping the gun and analyzer re- 
gions. If the system is not differentially pumped, careful attention must be given to 
baffling and shielding. In this mirror machine, the plasma flows out the ends of the 
magnetic bottle, and a weak background plasma may exist near the gun and the analyzer. 
This problem may be unique to this experiment because the magnets and the Penning 
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discharge are contained in a single large vacuum tank. 

In this experiment, the ion gun was run without differential pumping. The range of 

the accelerating voltage on the gun was 6 to 40 kilovolts. The background neutral gas 

9 12 

density in the discharge was 4x10 to 1. 7x10 neutral particles per cubic centimeter. 

An optical baffle was placed in the path of the beam to the plasma with a small opening 

for the primary beam. This eliminated radiation as an ionizing mechanism in the ion- 

gun bell jar. However, even with an optical baffle and metal screening in the pumping 

path between the ion gun and the vacuum tank, a weak plasma continued to diffuse into 

the gun bell jar. This usually occurred at background neutral gas densities of 8. 6x10 ^ 

neutral particles per cubic centimeter and anode voltages above 10 kilovolts. This weak 

plasma caused continuous arc discharges in the region of the gun. These discharges 

were reduced in intensity or disappeared as the discharge pressure, the anode voltage, 

or the gun accelerating voltage was reduced. 

The level of secondary signal strength in the plasma of interest should also be 

9 

examined. Estimates of the maximum density obtainable in this plasma are 10 to 

10 11 electrons per cubic centimeter at maximum density conditions. These densities 

are less than those seen in previous ion beam probe work, in which nanoamperes of 

secondary current were detected for microamperes of primary current in plasmas 

11 13 

where the densities were 10 to 10 electrons per cubic centimeter (refs. 1, 15, 
and 21). Here, we expect to see secondary currents 1 to 100 times smaller than those 
seen in previous work. Such currents were detected but only by nonautomatic operation 
of the probe. 

In an effort to understand more about the potential profile in the modified Penning 

discharge, we used an actuated Langmuir probe assembly to examine the region along 

a radial line through the midplane between the electrode rings, The results are shown 

in figure 14. The magnitudes of the measured profiles are neither those expected from 

the vacuum field solution nor those obtained from the computer modeling procedure. 

Results were obtained to within 3 millimeters of the anode, and the measured probe 

voltage was never greater than 0. 4 V„ and varied with pressure. The probe could not 

a 9 

measure the potential at the lowest discharge pressure, ;less than 4x10 electrons per 
cubic centimeter. Primary-beam orbits calculated from these potentials in the com- 
puter program were not in agreement with those measured. The ion beam results 
are considered to be more reliable than the Langmuir probe results, since the latter 
apparently disturbed the plasma to such an extent that profiles consistent with other data 
could not be obtained. 

The experimental results from this device are valuable in that they verify a partic- 
ular model of the discharge. The interpretation of the results is model dependent, but 
we think that the general characteristics of the size and shape of the potential would 
result no matter what initial model was chosen, as long as the boundary conditions of 
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Radius, cm 


Figure 14 - Potential measured with Langmuir probe in anode ring 
midplane. 

the modified Penning discharge were taken into account. Secondary beams have been 
observed, and a more extensive model with more parameters can be constructed. Since 
the model would have to describe the primary and all the secondary trajectories that 
were observed, confidence in the final model results would be increased. 

CONCLUSIONS 

It is possible to do ion beam probing of a plasma, the floating potentials of which are 
comparable in magnitude to the energy of the probing ion beam, and to obtain radial 
time-averaged potential profiles. The presence of large floating potentials in the 
plasma makes it necessary to use an iterative computer program to interpret data pro- 
vided by the primary and secondary ion beams. Several examples were observed for 
which the computer-derived radial potential profile correctly predicted the energy and 
momentum of both the primary and secondary ion beams. The experience obtained in 
developing mechanical and electronic components to probe plasma potentials compar- 
able to the primary-beam energy has been documented. This experience may be useful 
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in developing advanced ion beam probing systems for fusion-related plasmas which have 
strong electric fields or which contain high floating potentials. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 1, 1976, 

506-25. 
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APPENDIX - SYMBOLS 


B 

Cp C2 
E 

E i 

E 0 

e 

*tP *t2 
l 

M 

m 

N e 

N 0 

P 

q 

R a 

R g 

R i 

R o 

R r 

r 

s 

V a 

V ea 

V 

max 

v b 

0 ci 

0 mi 

<av)/v b 


magnetic field 

constants in assumed potential profile (eq. (6)) 
electric field 

ion energy along beam path 
initial ion energy 
electronic charge 

primary and secondary beam currents 

resolution length of electrostatic analyzer 

minimization parameter (eq. (7)) 

mass of particle 

electron density 

background neutral gas density 

pressure 

charge of particle 

anode radius 

radius of grounded boundary 
inner radius of potential maximum 
outer radius of potential maximum 
outer radius of potential minimum 
radius 

length along beam path 
anode voltage 
voltage across analyzer 
maximum potential 
velocity of beam ions 
calculated trajectory angle 
estimated trajectory angle 
effective cross section 



cp potential 

<P in minimum potential along plasma radius 

<P max maximum potential along plasma radius 

<p s space potential along beam path 
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